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Abstract 

Recently it has been claimed that ordinary perturbation theory (OPT) gives incorrect weak coupling expansions for 
lattice O(N) non-linear sigma models in the infinite volume limit, and in particular that the two-dimensional non- 
abelian models are not asymptotically free, contrary to previous findings. Here it is argued that the problem occurs 
only for one-dimensional infinite lattices, and that in general, OPT gives correct expansions if physical quantities are 
first computed on a finite lattice, and the infinite volume limit is taken at the end. In one dimension the expansion is 
sensitive to boundary conditions because of the severe infrared behavior, but this is not expected to happen in higher 
dimensions. It is concluded that spin configurations which are far from the perturbative vacuum have too small a 
measure in the path integral to invalidate OPT, even though they are energetically allowed for non-zero values of the 
coupling. 



Two-dimensional spin systems have been important 
field theoretic laboratories because of their similarities 
to four-dimensional gauge theories. For example, when 
N > 2 the O(N) symmetry is non-abelian and the the- 
ory has been shown to be asymptotically free [1,2]. The 
sigma model could also be regarded as a testing ground 
for perturbation theory. As for gauge theories, the per- 
turbative vacuum, with all spins aligned, is quite dif- 
ferent from the true vacuum, since there is no magne- 
tization for d < 2 [3]. A further similarity is that the 
expansion proceeds in powers of the coupling and the 
fields, although there is no mass to damp out configura- 
tions far away from the perturbative vacuum. Therefore 
one might worry that very long wavelength excitations, 
where the fields eventually get large, will contribute sig- 
nificantly to the path integral, and that ordinary pertur- 
bation theory (OPT) may not account for them correctly. 

Thusly Patrascioiu has argued that OPT gives incor- 
rect low-temperature expansions for the free energy and 
spin correlations of the O(N) non-linear sigma model, 
for N > 2 and infinite volume lattices [4]. (For finite 
lattices there always exist temperatures low enough so 
that the spins are all relatively aligned, and OPT should 
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have no problem.) To remedy this, he formulates a new 
perturbation expansion which treats the gradients of the 
fields, rather than the fields themselves, as small quan- 
tities. Despite the method's computational difficulty (it 
is very nonlocal), he extracts a value for the one-oop 
spin-spin correlation function in two dimensions, which 
yields an opposite sign for the Callan-Symanzik /3 func- 
tion, relative to the OPT result. 

It is important to investigate this claim, for if it is true, 
then it could conceivably be that the /3 function for QCD 
or other gauge theories is different — in magnitude, if not 
in sign — from what is presently accepted. In this letter I 
show that OPT gives the correct result in one dimension, 
where the theory is exactly soluble, ;/ it is formulated 
on a finite lattice with the correct boundary condition. 
Since this result is independent of the lattice size, L, it 
is trivially correct in the infinite volume limit. This pro- 
cedure should give correct results in higher dimensions, 
where the IR behavior is milder. It is shown that starting 
with an infinite lattice gives results consistent with the 
exact solution for N — > oo (as well as the infinite volume 
limit of the finite lattice calculations) except in the one- 
dimensional models, apparently because of their more 
severe IR divergences. 

The controversy over OPT's validity stems from the 
difference between the OPT result for the free energy 
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density on a <f -dimensional, infinite lattice [5], 

^^ + I(Ar-l)ln, + ^ + 0(^) (1) 
versus the expansion of the exact result for d = 1 [6] 
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where the lattice partition function is given by 



r ( ] d ) 

Z= f[dS5(S2-l)exp -J]2 S '' S ^ 

^ X \8 X (1=1 



(3) 



OPT proceeds by rewriting S as (tt, <x), where n is 
(N - l)-dimensional, and solving the constraint, so that 
cr = ±(1 - n 2 ) 1/2 . The vacuum state is chosen to be 
cr = 1, and configurations with cr < are discarded, as 
they give, naively, contributions of 0(exp(-l/g)). After 
rescaling n 2 to gn 2 , Z becomes 



Z = exp(dL d /g) 
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exp 



(4) 



To do perturbation theory, the integration region must 
be extended from n 2 < 1 /g to n 2 < oo and the radicals 
Taylor-expanded. 

More crucially, the zero modes of the n field must be 
removed. This could be accomplished by introducing a 
magnetic field in the direction (i.e., n = and cr = 1), 
which is equivalent to a mass for the n, and removing it 
at the end of the calculation. However, this method is 
known to give wrong results even for a two-spin sys- 
tem, beyond the tree level [7]. It is easy to see why 
this happens, heuristically: when the factor (1 - gn 2 ) 1 ^ 2 
in the measure is reexpressed as a term in the action, it 
contributes a mass of the wrong sign, m 2 = -g. Then, 
as the magnetic field becomes smaller than g, the ac- 
tion becomes unbounded from below, and the gaussian 
integrations no longer make sense. A better way is to 
use the global O(N) invariance to fix one of the spins, 
say So, to be in the direction. More generally, one 
can use the Faddeev-Popov procedure to fix any linear 
combination of the S x , thus removing the zero modes, 
and this gives agreement with exact solutions in every 
case that has been checked [7]. However, fixing spins is 
not sufficient for an infinite lattice, because it provides 
a low-momentum cutoff of 0(1 /L) (L being the lattice 



size) which vanishes as L — > oo causing the n propaga- 
tors to become undefined. Thus one is forced to bring in 
a magnetic field again. Since this is a bad procedure for 
finite systems, there is no a priori reason for it to work 
on an infinite lattice. But this is exactly how (1) was 
obtained; hence the discrepancy with (2). 

Patrascioiu has a different explanation for why OPT 
goes wrong as L — > oo in 1 dimension, however. Note 
that the two-point function is given by 



(S • S x ) = g{n ■ n x ) 

+<(l-«ag) 1/2 (l-«^) 1/2 > 
= l+g(N-l)D(x) + 0(g 2 ) 



(5) 



at tree level, where D(x) = G(x) - G(0) and G(x) is 
the massless scalar propagator in d dimensions. The 
large-x dependence of D(x) is -\x\, — In |x| and 
for d = 1,2,3, respectively. Therefore we have con- 
figurations in which the spins wander far away from 
en, the perturbative vacuum, over a distance \x\ ~ g~ l 
(\x\ ~ exp(^ -1 )) for d = 1 (d = 2), whereas these are 
energetically suppressed for d > 3 (magnetization). Of 
course, OPT is not designed to account for spins point- 
ing far away from e^, since cr < was discounted. The 
argument is that for finite L there should be no problem, 
because there always exists a g small enough so that gL 
or g In L is «: 1 : then no large excursions away from 
ejv are allowed. For L = oo, d < 2, the phase transition 
occurs exactly at g — 0, so that for finite g these long 
spin waves are unsuppressed, and OPT may fail. 

In order to clarify why OPT fails in one dimension, 
I compute the correlation function (So ■ S x ) in d dimen- 
sions, using a somewhat different method, which is to 
obtain the non-linear cr model from the infinite mass 
limit of the linear model. This has the advantage of 
avoiding the OPT approximations which made the long 
spin waves impossible to represent correctly. It is in- 
structive to do the calculation on both finite and infinite 
lattices, for it will become apparent that only in one di- 
mension does any discrepancy appear. 

First consider the finite lattice. We start by noticing 
that the 6 functions in (3) can be rewritten using the 
identity 6(x) = lim M ^oo(M/ V^) exp(-M 2 x 2 ). After the 
change of variables S — > ( s/gn, 1 + yfgcr) and "gauge 
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Figure 1 : Propagators and vertices for the linear cr model, including the lowest order vertex from the Faddeev-Popov determinant 



fixing," the partition function is 
exp 
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where 6(F) removes the zero modes of n, and A F 
is the associated determinant. A convenient choice 
which preserves the n propagator's translational invari- 
ance is F = Yjx^x, in which case A F turns out to be 
(1 + ^L- d Z x o- x ) N - 1 (seeref. [7]), and the bare prop- 
agator is 



{Jtx ■ no) 



(0) 



(N-l)G(x) 

d v exp(i2np-x/L) 

" " ' ~ P(P) ' 



pi=Q 



Pip) = 4 Yj s in2 "P»I L - 



(7) 



The propagators and vertices of this theory are shown 
in fig. 1, including the lowest order vertex due to A F . To 
compute the two-point function 

<S • S x ) = 1 + yjg((a x ) + (ao)) 

+ 8({o-qo- x ) + {n Q ■ n x )) (8) 

to 0(g 2 ), we need the diagrams listed in fig. 2, whose 
leading M 2 -dependences are nonvanishing. Fig. 3 gives 
examples of graphs which vanish as M — > oo because 
they have more cr propagators than non- A F vertices. Al- 
though some of the graphs in fig. 1 diverge like M 2 , their 



sums are finite, as expected, and the factor M Ld in Z is 
exactly cancelled by [det(-A 2 +M 2 )]~' coming from the 
cr integrations. After much algebra, the one-loop con- 
tribution to (So ■ S x ) obtained is 



(S -S X )W = g 2 [-(N-l)(N-2) 



(9) 
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in exact agreement with the OPT result of Hasenfratz 
[7]. This just shows that the OPT approximations — 
ignoring cr < configurations, and enlarging the in- 
tegration region beyond the radius of convergence of 
(1 - g7T 2 ) 1/2 's Taylor series — are justified for L < oo 
as expected. 

Before carrying out the finite lattice calculation of 
(So • S x ), I will show that the procedure leading to (9) is 
correct for d = 1 and any L (and thus, as L — > oo). This 
case is important because: (1) only here is the exact so- 
lution known (unless N = oo); (2) the low-energy spin 
waves should affect OPT most severely when d = 1, by 
Patrascioiu's argument; and (3) the only known failure 
of OPT occurs when d = 1 . For simplicity, let us only 
consider the average interaction energy per spin, 
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From (9) we find that 

= i-i^AT-i) 

+ ±g 2 (N-l)(N-5) + 0(g 3 ) (11) 
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Figure 2: Diagrams contributing to (So • S x ) to second order, with leading M dependence greater than M~ 2 
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whereas the expansion of the exact result is [6] 

/exact = 1-±S(W-1) 

+ ig 2 (JV-l)(jV-3) + 0(g 3 ) (12) 

independent of L. In comparing these it must be re- 
membered that the exact result is obtained for an open 
chain of spins, in which the ends do not interact with 
each other, whereas OPT is normally done on a periodic 
lattice, which is the case in (1 1). To repeat the OPT cal- 
culation (4) for the open chain, one must use the open 
chain lattice propagator, which in Yj X n x = gauge and 
at 0(g°) is exactly 

G open (x,y) = G Xi y + L(G xL - G x \){G yL - G y \). (13) 

Here G xy is the periodic lattice propagator in (7), which 
for d = 1 has the closed form 

G^ = ±{L-L- l )-\\x-ymo&L\ 

+ ^(i-ymodi) 2 (14) 

Therefore the momentum sums in loop graphs can be 
done explicitly. The result is (12), including the ex- 
act L-independence which distinguished it. Thus, OPT 
agrees with the exact result in d = 1 as long as the same 
boundary conditions as in the exact solution are used. 
This does not mean that (1 1) is incorrect. To prove that 
the difference between (11) and (12) is due to the phys- 
ical difference made by the extra link interaction, and 
not some fluke of OPT, I applied Patrascioiu's method 
to the periodic spin chain with L spins and confirmed 
(11). (This letter has no complaint with the method of 
ref. [4], but only its result.) 

Of course, the exact expressions for / must agree for 
the two boundary conditions as L — > oo. The only plau- 
sible explanation for the difference in the perturbative 



expansions is that the exact expressions differ by a func- 
tion like 

/open - /elosed ~ L~ X SX P [-g Z L(N - 1)(N - 2)/12] (15) 

for example. Even though this vanishes as L — > oo, its 
expansion in powers of g does not, and in fact it pre- 
dicts that the weak-coupling expansion of / c i ose d is IR- 
divergent, starting in 0(g A ). From the point of view of 
the perturbative calculation, the reason for sensitivity to 
the boundary condition, which is a 1/L effect, is that 
the propagator is linear in distance and so acquires val- 
ues of O(L) in the loop diagrams. In two dimensions 
the propagator is logarithmic, so we might expect the 
L in the exponent of (15) to be replaced by InL. The 
resulting expression cannot lead to differences in per- 
turbation expansions, so the effect explained by (15) is 
almost certainly unique to d — 1 . 

Having seen that OPT on a finite lattice gives the cor- 
rect expansion for the energy density as L — > oo, in one 
dimension, we now examine what happens when L is 
taken to be infinite at the outset. In this case, as noted 
previously, the Faddeev-Popov terms in (6) are insuf- 
ficient for regulating the n propagator; the field must 
be given a mass, //, to be eventually removed. But this 
spoils the O(N) symmetry, and so one is not justified to 
insert gauge-fixing terms: A F 6(F) should no longer ap- 




Figure 3: Examples of graphs which vanish asA/-> oo. 
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pear in (6). (Actually, for the gauge choice F = Yix n x 
it turns out not to matter whether one keeps the 6(F) 
factor, when L = oo.) It is straightforward to show that 
the only effect of these changes on the result, (9), are 
(1) to replace the second term by its L — > oo limit (re- 
call that D(x) is IR-finite); and (2) to replace the first 
term, the A F contribution, by -^g 2 (N - l)\x\ 2 (zero) 
when d = 1 (d > 1). This agrees fully with the OPT 
calculation on an infinite lattice for d — 2 which was 
done by Elitzur [2], Since the present method was de- 
signed to correctly measure the contributions of the long 
spin waves, yet it agrees with OPT for d > 2, we con- 
clude that the long spin waves do not cause OPT to give 
incorrect results when d > 2. Furthermore, the problem 
in d = 1 is clearly seen to be due to the noncommutativ- 
ity of the two limits g — » and p — > (where by g — > 
I mean developing the asymptotic expansion), just as 
was suggested earlier. On the other hand, this analysis 
shows that these limits do commute for d > 2. (More 
precisely, the order p — > 0, g — > 0, L — > oo commutes 
with L — » oo, g — »0, p — > for d > 2.) 

This conclusion derives further support from compar- 
ison with the N - oo limit of the model, where again 
exact solutions are known [8]. This limit exists if the 
coupling is rescaled so that f3 = 1 /gN is held fixed as 
N — > oo. Then, for example, (12) would become 



d=l Ip op 



(16) 



and using the fact that D(\) = -± for d = 2, (5) and (9) 
give us 
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There should be no problem with first doing the small g 
expansion and then taking the L — > oo limit, since the 
latter requires g to be infinitesimal. Indeed, the exact 
solution for = oo, in one dimension [6] is 



f N - 



(18) 



d=i 1 +(1 +4/3 2 ) 1 / 2 

which agrees with (16) when expanded in powers of 
1 1/3. In two dimensions, the exact solution [9] is 

i-X 
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where m is a dynamically generated mass, which for 
large /3 is given by 

m 2 = 32 exp(-47r/?) + 0(exp(-8^8)) (20) 



For large (3, hence small m, the integral is dominated by 
its p = contributions, and we can approximate it by 
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= j6 - — ln(32/m 2 ) + 0(exp(-47r y 6)), (21) 
An 

where A = (9(1). The important thing to notice is that 
when this is substituted into (19), the expansion of / 
in powers of 1//3 terminates after the 1//3 term. This 
agrees with the OPT result, (17), whose 1//3 2 term is 
seen to vanish. In contrast, Patrascioiu's result for the 
energy density in d = 2 contains a piece of order g 2 N 2 
which when reexpressed in the large-Af limit is a 1 10 2 
term. There seems to be no subtlety of ordering limits 
in this comparison, since (19) was derived for an infinite 
lattice, and as has already been mentioned, g — > and 
N — > oo are expected to commute. 

Finally, the /3 function calculation of Polyakov [1] 
merits attention because it is on a completely indepen- 
dent footing from the OPT method, yet it gets the same 
answer, /3(g) = -g 2 (N - 2)/2n + 0(g 2 ). The salient 
point is that it uses Wilson's method of integrating out 
the high-momentum components of the spin field and 
seeing how the coupling is renormalized in the effective 
action for the remaining low-momentum components. 
Thus, /3(g) is determined by the short wavelength fluctu- 
ations, and is insensitive to the long spin wave effects by 
which Patrascioiu explains the difference between his j3 
function and the standard result. 

In conclusion, I have shown that OPT on a finite lat- 
tice is capable of giving the correct asymptotic weak- 
coupling expansion for d = 1, even in the infinite vol- 
ume limt; therefore the same procedure should work 
in two dimensions, corroborating the standard result: 
asymptotic freedom for N > 3. The fact that the 
M — > oo limit of the linear cr model (performed here), 
the known N - oo solutions, and Wilson's renormaliza- 
tion group (ref. [1]) lead to the same answer gives one 
yet more confidence in ordinary perturbation theory, as 
applied to the O(N) nonlinear cr models. 

I am indebted to John Preskill for suggesting this 
problem and for many stimulating discussions, and also 
to J. Feng and B. Warr for the latter. Lately I learned that 
H.D. Politzer and G. Siopsis were carrying out related 
investigations [10]. I thank them for helpful exchanges 
of ideas. 
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